Tilt Texture Domains on a Membrane and Chirality induced Budding 
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We study the equilibrium conformations of a lipid domain on a planar fluid membrane where the 
domain is decorated by a vector field representing the tilt of the stiff fatty acid chains of the lipid 
molecules, while the surrounding membrane is fluid and structureless. The inclusion of chirality 
in the bulk of the domain induces a novel budding of the membrane, which preempts the budding 
induced by a decrease in interfacial tension. 
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The budding of membrane vesicles and subsequent fis- 
sion is the only means of traffic between compartmental 
boundaries within a living cell Budding of vesicles 
clearly requires membrane deformation — coat proteins 
such as Clathrin and COPs I-II have been postulated to 
play a crucial role in facilitating this deformation. How- 
ever the precise role of these coat proteins in budding 
events remains poorly understood 

Membrane deformation leading to budding may also 
be triggered by changes in the structure and composi- 
tion of lipids. An oft-quoted mechanism, especially in 
the context of artificial membranes, is the phase segrega- 
tion in multi-component systems which leads to budding 
due to a decrease in the interfacial energy ||,[| between 
the lipid components. However the resulting bud has a 
size Rbud > 2k/ A, where k is the bending modulus of the 
membrane and A the interfacial tension, which for typ- 
ical lipids is of order 5/im, much too large compared to 
the budding events in the cell. In this letter we propose 
a novel budding mechanism in a multi-component lipid 
membrane, where one of the species is chiral and acquires 
a tilt. This chirality-induced budding can result in a bud 
size as small as 500A. 

Consider for simplicity a lipid bilayer membrane, where 
the upper layer is built of two species, a majority compo- 
nent (e.g. DPPE) and minority component (e.g. sph- 
ingolipid), while the lower layer consists of a uniform 
density of DPPE Let the minority species (sphin- 
golipid) acquire a tilt with respect to the local normal of 
the membrane. The tilt || could be intrinsic (when the 
temperature T < T m , the main transition temperature 
of the chosen lipid) or acquired due to association with a 
third component (e.g. cholesterol). We will work in the 
temperature regime where there is at least a microphase 
separation of the two lipid components over a length scale 
R of the order of say 500A. Note that most cellular lipids 
are chiral and so, on acquiring a tilt, this molecular chi- 
rality will express itself at length scales much larger than 
molecular dimensions. The chiral strength may of course 
be enhanced by forming complexes with cholesterol . 

To start with, let us assume that the sphingolipid com- 
plex is spread, with a fixed, uniform density, across a 
domain of area A, small enough so that over this scale 



the membrane may be assumed flat. As a result, our 
description entails a decoration of a planar domain by 
a two-dimensional (2d) unit vector field m, representing 
the projection of the stiff lipid chains onto the plane of 
the membrane [p|, while the surrounding membrane is 
fluid and structureless. 

Our low energy hamiltonian for the membrane tilt 
domain is given, to quadratic order in the fields by 
E = Eb + El ||, where the bulk energy 

E B = I fci(V • m) 2 + fc 2 (V x m) 2 + fc c (V • m)(V x m), 

J A 

(1) 

and the interfacial energy, 



Er 



dl (ao + (Ti(m • n) + er 2 (m x n)) 



(2) 



The unit normal to the boundary n aims everywhere into 
the domain. Note that V x m is a pseudoscalar in 2-dim 
and fc c (V • m)(V x m) is the bulk chiral contribution to 
the energy. The magnitude of k c is set by the density of 
cholesterol present in the sphingolipid complex, ao, the 
isotropic line tension, is of the order of 10~ 8 dyne. The 



Frank constants k\ and ki are ~ lOfc^T,., 



while the 



interfacial energy of a domain of radius 0.05/im is of the 
same order ]l0|]. This implies that thermal fluctuations 
can be ignored at the temperatures of interest. We set 
k\ = 1 as our unit of energy. 

The hamiltonian (|l|), (0) has been studied in HfyJ in 
the context of textures on Sm-C* films and Langmuir 
monolayers respectively. Apart from the context, there 
are three major points of departure — (i) our variational 
calculation and Monte Carlo (MC) simulations for the 
T = phase diagram are 'exact' giving results very differ- 
ent from |llj , (ii) we study the effect of strong anisotropic 
line tension and bulk chirality on the shape and texture 
of the domain and (iii) we couple the tilt field to the lo- 
cal curvature of the membrane and study the effect of 
chirality on membrane shape. 

We rotate our coordinate axes and redefine the coef- 
ficient of anisotropic line tension to obtain the simpler 
expression |,|| E L = § dl (a + cr(m • n)). The T = 
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phase diagram is obtained by minimising E variationally 
subject to the constant area A constraint. In all cases, 
we have corroborated the results of our variational cal- 
culation by an MC simulation using simulated annealing 
to avoid getting stuck in local minima. In this sense our 
variational calculation is 'exact'. 

We first study the phase diagram when the isotropic 
line tension is extremely high forcing the domain bound- 
ary to be circular. To begin with, let us turn off bulk 
chirality (k c = 0) and set fc 2 = k\. The T = phase dia- 
gram in the R — a plane (where R = ^ A/ir) has the fol- 
lowing 4 phases separated by first-order boundaries (Fig. 
1) : (1) Boojum phase where the centre of the circular 
domain is a distance a away from the core of a virtual 
defect of strength N. Placing the origin of polar coordi- 
nates (r, 9) at the core of the boojum, we can describe 
the texture m = (cos <fi, sin0) by the equation 

</> = N9 + ar sin 9 + c 2 r 2 sin 29 + c 3 r 3 sin 30 + . . . . (3) 

The energy E is minimised with respect to the variational 
parameters (N, a, {c„}) for fixed (R, a). We find that the 
optimal value of N is exactly 2, independent of a. In 
drawing the phase diagram we parametrise the texture 
by N and c\ alone; inclusion of c 2 and terms of higher 
order lowers the energy by 1% at most. Fixing the area 
at say R = 1 (in units of a microscopic cutoff length r c ), 
an increase in a leads to an increase in c\ and pushes the 
core of the virtual defect towards the domain centre, till 
at a — 1.7 we pass into the annular phase. 
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FIG.l. Phases of the tilt texture domain with circular 
periphery: (1) Boojum (2) Achiral Annular (3) Chiral 
Annular (4) Hedgehog with e c = 0. Arrows indicate di- 
rection of m field. 

(2) Achiral Annular phase with inner and outer radii 
ri and r 2 respectively. At the outer rim of the an- 
nulus m is directed radially outward while at the in- 
ner boundary m is inclined at an angle a to the local 



normal n. With the origin of polar coordinates at the 
centre of the annulus, the texture may be described by 
4> = 9 — a(r 2 — r)/{r 2 — r\ ) , where r\ , a are variational pa- 
rameters. The inner radius r\ monotonically decreases as 
we increase a such that m always points radially outward 
(a = 0) at every point of the domain. As one crosses the 
first order phase boundary (Fig. 1) we encounter a new 
annular phase with a chiral texture. (3) Chiral Annular 
phase where a jumps to an angle greater than 7r/2 at 
both boundaries. Consequently the bulk texture is chiral 
with the order parameter C = J d 2 x (V • m)(V x m) ^ 
even though there is no bulk chiral interaction. This chi- 
ral phase is doubly degenerate and can spontaneously 
acquire either sign. (4) Hedgehog phase with a defect 
at the domain centre and a texture described by (f> = 9. 
The energy includes a core contribution e c of the defect 
of size r c . 

The transitions between phases 1 — > 2 and 2^3 arc 
discontinuous (indicated by the change in slope of the 
energy branches as a function of a) and weaken as the 
domain size R shrinks. 

As cr is reduced, the boundary of the domain is 
no longer constrained to be circular. We may then 
parametrise the boundary by smooth deformations of a 
circle, r(9) where the origin is at the centre of the circle 
of radius ro and < 9 < 2ir is the angle to the polar axis, 

r{9) =r Q (^l + a n cos(n0)^ . (4) 

The amplitudes {a n } are additional parameters in our 
variational calculation. We find that the texture through- 
out the R — o~ plane is a boojum with the domain bulging 
out near the equator and flattening near the poles as a 
increases from zero (henceforth we shall set cr = 1> as a 
result our unit of length becomes of order 0.1 ^m). These 
shapes and textures win over the annular and hedgehog 
phases. 

We now turn to a description of phases with nonzero 
bulk chirality k c , which without loss of generality we take 
to be positive. In the phase diagram, Fig. 2, we have 
set cr = 0: (1) Uniform phase, when \k c \ < 2, where 
the domain is circular and the tilt field m is uniform. 
(2) Spiral Hedgehog phase where the domain is circular 
with a texture such that m at every point is inclined at 
an angle a to the radially outward ray emanating from 
the centre of the domain. The radial and the tangential 
components of m in this spiral phase are m r = cos a, 
me = sin a, while the energy is given by 

E = 2n(r c + R) - sin 2a - 1) log — + e c . (5) 

Z r c 

The optimum value of a is 7r/4 when \k c \ > 2. This phase 
clearly gives a nonzero value for the chiral order param- 
eter C. (3) Tweed phase domain is again circular though 
the texture loses the circular symmetry of the previous 
phase. Our ansatz for the texture is motivated by the re- 
sults of an MC simulation of 3055 particles carrying 0(2) 
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spins on a triangular lattice with a hamiltonian obtained 
by discretising E, for a — 0, k c = 25.5. We have per- 
formed a simulated annealing from (UbT) -1 = 0.1 — > 300 
starting from a variety of initial conditions |Q| to obtain 
textures as in Fig. 2 (lower inset). 




FIG. 2. Phases of the tilt texture domain with bulk 
chirality and a = : (1) Uniform (2) Spiral Hedgehog 
(e c = 0) (3) Tweed. Top inset shows variation of the en- 
ergy branches as a function of k c at R = 1 indicating dis- 
continuous transitions. Notice that the energy decreases 
rapidly in phase (3) as k c increases. Lower inset shows 
the texture within a unit cell ; this pattern repeats peri- 
odically to form the Tweed phase (3). 

A variational ansatz may be easily written down for such 
a texture: place the origin of coordinates on the bound- 
ary of the domain and parametrise the texture by a 
length, I, such that m x = | cosa;/Z| and m y — — | sinx/i). 
This texture confers a net chiral strength to the domain 
— (V ■ m)(V x m) has the same sign everywhere, though 
the sign of individual terms, V • m and V x m, varies from 
one stripe to the next. The energy of the tweed domain 
is 

k r 2R 2x 1 

E = 2itR--^I dx\J2xR~x 2 | sin — | + -j . (6) 
21 Jq it 

The variational parameter I approaches zero to minimise 
the energy, corresponding to infinitely many tweeds. 
Higher order derivative terms in the hamiltonian would 
however cutoff this monotonic decrease at some scale I* 
& 

On turning on the anisotropic line tension er, the 
boundary of the domain deviates from a circle giving rise 
to a variety of new phases. We will make a detailed study 
of the phase diagram in a later paper. Let us move on, 
for the moment, and allow the membrane containing the 
domain to be flexible, thus promoting a coupling between 
the tilt field m and the curvature tensor Kij Q . At the 



same time we have to treat the membrane as a bilayer: 
we thus define the curvature, tilt and lipid density fields 
on both the upper and lower leaves of the bilayer and 
project these variables onto the neutral surface of the 
membrane 0]. In covariant form, the bulk energy to 
lowest order in derivatives, is given by [fl5f , 

Eb = I [fci(divm) 2 + fc2(curlm) 2 + fc c (divm)(curlm) 

J A 

+c Q Kl + k (Kj) 2 + (3m l m 1 K) + c*j l:j m k 'm l k{\ , (7) 

where the divergence and curl are written in terms of 
the covariant derivative Di as, div m = Dim 1 and 
curlm = '-ffDj-m 1 . The appearance of 7.^ , the com- 
pletely antisymmetric tensor on a 2d curved manifold 
|13| , indicates that the coupling of the texture to the 
membrane curvature also depends on the chiral strength 
of the sphingolipid complex (?| in the domain. To high- 
light the effect of chirality, we set Co = (3 = and only 
include the isotropic line tension in the interfacial energy 
El 0. 

We will see that the effect of increasing Cq for fixed size 
R is to induce the domain to leave the plane and form 
a spherical bud of radius R/2 connected to the plane by 
an infinitesimal neck ! More significantly, at fixed but 
large enough Cq this chirality induced budding occurs at 
a much smaller size R than the phase separation induced 
budding caused by a decrease in the interfacial tension 
between the two lipid components 

For instance if we fix k c at phase (2) of Fig. (2), we 
find that the membrane remains planar and the texture a 
spiral hedgehog for small Cq (Fig. 3). A further increase 
in Cq leads to a discontinuous budding transition where 
the texture decorates a spherical bud with a vanishingly 
thin neck. Both the chiral terms, proportional to k c and 
Cq, gain in energy by this budding, even at the cost of 
producing a defect at the pole. This may be seen by 
setting the latitudinal and longitudinal components of m 
as mg — \j\[2 and — 1/^/2, a choice which smoothly 
goes over to the spiral hedgehog phase in the limit of the 
planar membrane. The membrane shape is parametrised 
as a sphere of radius r attached to the plane by a neck of 
radius rQ. The energy of this conformation is El = 2itro 
and 



E r 



R 2 



R 2 



ir(k c 



*~ 6 ° cos 2 

de c ^l (8) 

smO 



in addition to a defect core energy e c , where 8 C = r c /r 
and ro = rsin^o- We have ignored the small curvature 
energy cost at the junction where the spherical bud meets 
the plane. Variation suggests that the energy is min- 
imised when r — > and r = R/2 (Fig. 3). 

At larger values of Cg , the spherical bud gives way dis- 
continuously to a cylindrical tubule (Fig. 3). We see that 
the energy is minimized when the domain wraps around 
a cylinder of radius r and length L attached to the plane 
at one end and capped by a hemisphere with a polar 
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defect core at the other (Fig. 3). The m field on the 
spherical cap is as described earlier; the axial and tan- 
gential components of m on the cylinder are m z = 1/V2 
and = l/v2- This choice corresponds to a texture 
where the lines of m wrap around the cylinder as a helix 
of pitch proportional to r Wfa. 




FIG. 3. Phases of the tilt domain coupled to membrane 
shape with k c = 2.5, a = and k — 5 : (1) Planar (2) 
Spherical Bud (3) Cylindrical Tubule. With these pa- 
rameters, domain induced budding due to a gain in line 
tension j^] occurs at R c = 10. Inset shows the variation 
of the ratio a = L/r with Cq for a domain of size R = 1. 

Note that m varies smoothly across the junction of the 
cylinder and the hemisphere, and extrapolates smoothly 
to the spiral hedgehog texture in the limit of the planar 
membrane. The energy of this conformation is El = 2ttv 
and 



Er 



2tt 



(k - c* r) + 7T (| - c o r ) 
cos 2 9 



2r 2 



1 



f7T — Oq 

-7t(> c -2) / dO 



sin ( 



(9) 



in addition to a core energy e c (we have again ignored the 
curvature contribution coming from the neck). The op- 
timal L/r obtained variationally, increases rapidly with 
c* (Fig. 3 (inset)). 

As seen from the inset of Fig. 2, the energy plummets 
as k c increases. This implies that if k c was chosen such 
that phase (3) of Fig. 2 were the equilibrium texture at 
Cq = 0, then the domain would gain significantly in bulk 
energy by budding, first into a spherical bud and then 
into a tubule. The range of Cq for which the spherical 
bud obtains would however be reduced. 

The variational calculation for the T = energy- 
minimising shapes of a flexible, fluid membrane having a 
tilt domain is of course not exact. A more general varia- 



tional calculation will at most shift the phase boundaries 
slightly — in any case the prediction of a chirality-induced 
budding at small R is robust and should be testable in ex- 
periments on artificial membranes consisting of carefully 
chosen lipid components. More importantly, we believe 
that this may have some relevance in the budding of cell 
membranes which are not assisted by coat-proteins, as in 
the endocytosis of raft components |18|. 
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